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ABSTRACT 


A fractal tree of nested viscoelastic boxes is proposed to describe the elastic after-effect in some 
composite solids that show evidence of a hierarchical structure, such as clay, some rocks, 
concrete, and even some biological tissues (as well as materials used in artificial grafts). A 
functional equation is derived and solved for the Laplace’s transform of the creep kernel of the 
whole hierarchy. Inverting the transform, Hopkinson’s power law is obtained as a fairly 
accurate approximation, after an initial time interval. The exponent in the power law is 
obtained as a well-defined function of the branching numbers and scaling parameters of the 
viscoelastic hierarchy. Then we consider a composite solid with both elastic (impact) strain 
response and power law type elastic after-effect, for an applied stress step. The stress response 
to an applied strain is given in terms of Rabotnov’s fractional exponential functions. The 
stretched exponential stress relaxation to an applied strain step is derived as a fairly close 
approximation. According to the present mathematical model, for a given composite solid, the 
stretch parameter of the stretched exponential should be equal to the exponent of the power law. 
The experimental lower bound (0) and upper bound (1) for this parameter are derived from the 
theory. It is shown how, at least in principle, these results could be used to design materials 
with given values of the power law and stretched exponential parameter. Finally, a generalized 
transmission line approach (a so called fractance approach) is developed to relate the strain 
and stress responses of hierarchically structured composite solids. A comparison is made 
between the results obtained using the fractance approach and the results obtained using the 
functional equation approach. 


Key Words: Solid composites, elastic after-effect, fractal hierarchy, power law strain, 
stretched exponential stress relaxation, fractance, functional equation, Rabotnov’s fractional 
exponential function. 


INTRODUCTION 


It is known that certain composite solids, that show evidence of having a 
hierarchical structure, like clays, some rocks, concrete and biological materials, often 
present hyperbolic delayed strain responses (known as power law creep) to a step of 
stress, and exponentially stretched stress responses to a step of strain, at least very 
approximately, during long periods of time, after a suitable initial time interval. 

The power law creep was discovered by W. Weber in 1836, studying the delayed 
strain €,(t) response of the silk threads of galvanometers under a constant tensile 


loado,. When a constant stress of magnitude o, is applied at the origin of the time 


variablet, the strain response é(t) to an applied step of stress has a linear elastic 
component €, and a delayed component ¢, (t): é(t) =E+é, (t) 


If J, is the elastic compliance modulus, the elastic response is given by & =J,:o> 


P 
The delayed strain seems to follow a hyperbolic law e,(t)= K {2] -O, Here K isa 
(a 


positive parameter. This power law fits very accurately experimental data for times ¢ 
greater than the characteristic timer, . The parameter p is always between 0 and 1. So, 


P 
: : 2%. t 
the whole strain response to a step in stress is given by: é(t)= c ot K- 4) J 
T 


For a variable applied stress o(t) the delayed strain response of a linear viscoelastic 


power law material is given by ¢,(t)= [ Q(t -1')-o(t')- dt 
P 
If o(t)=0, for t>0 we obtain ¢,(t)=0, [ Q(t')-dt' Then [ Q(t')-dt'= {=| so 
ic 


(l-a)-K 


€ 


the viscoelastic kernel is Q(t) = O(r, ): (Zs) with a =1—p and Q(z.) = 


If o(t) =o, for 0<t<T), and zero otherwise, then the delayed strain response for t > T, 


l-a la 
is: eo [0lr-r bara, I: C 2 | But in this case, if 


em 


t > +oothe delayed strain response ¢,(t)>0. As, of course, the instantaneous elastic 


response disappears also, the material behaves as a solid. This behaviour justifies the 
usual denomination of “power law solid” applied to this kind of materials. 
Now, notice that the power law response can be writtene, (t)= J (t)-o, 


with J(t)=J(1):t?.The compliance J(t)verifies the following functional equation: 
I(t, -x)=M(t,):F(x) Here M(t,)=t? and F(x)=J(1)-x?. 

With ¢, fixed, x €[0,1] andt=1,-x €[0,t,], this functional equation characterizes a 
scale invariant response J(t) in a time interval[0,t,]: M(t,) gives the amplitude and 


r (x) gives the shape (always the same) of the response. Conversely, the only regular 
solutions of the functional equation of scale invariance are the power functions [1]. 


[om 


P 
The stretched exponential stress relaxation o,(t) = o(0)-exp | (| was 
E 


discovered by R. Kohlrausch in 1847 while studying the relaxation of shear stress under 
conditions of constant shear strain, also in the silk threads of galvanometers. For times f¢ 
greater thanz, , the stretched exponential fits to experimental data fairly accurately. The 
parameter p also is always between 0 and 1. 


After its discovery in 19" century, both hyperbolic and stretched exponential 
responses were often found studying the mechanical properties of many bodies, 


including glasses [2], polymers [3], rocks [4] [5], concrete members [6] and samples of 
biological tissues [7] [8]. 

There is another fact that is closely related with the hyperbolic response. It is the 
relative constancy of the quality factor corresponding to elastic wave dissipation in the 
Earth, for most earth materials and for wave periods comprised between 10° and 10° 
seconds: it lies between 100 and 1000 [4]. In 1976, in the 46" annual meeting of the 
Society of Exploration Geophysicists (USA), Strick presented a phenomenological 
mathematical model that explains this remarkable behaviour in consolidated 
sedimentary rocks using a continuous fraction stemming from an infinite network of 
springs and dashpots with increasing elastic and viscous parameters away from the 
point of stress application [9]. 

In this type of model, strain and stress responses result from the addition of partial 
strains or stresses that correspond to the degrees of freedom in an infinite hierarchical 
structure [10], [11]. In order a given degree of freedom be able to relax, the precedent 
must relax before. As consequence, the time scale of evolution of a given level is 
subordinated to the time scale of the precedent level: each level responds slower than 
the precedent level but faster than the next one. 

The main purpose of the present paper is to generalize Strick’s model, applying the 
elementary mechanical theory of composite materials to a suitably defined hierarchical 
structure of elastic elements of different sizes and shapes immersed in a viscous matrix’. 
Under the assumption that the dominant viscous dissipation occurs at the interfaces 
between the elastic elements and the viscous matrix, we obtain the power law creep that 
describes the delayed strain response of the above mentioned composites. 

In order to do that, we derive a functional equation in Laplace domain, whose 
solution in time domain corresponds to power law creep. This derivation follows an 
approach closely related with Liu’s explanation of constant phase angle impedances and 
1/f noises in rough electrode interfaces [12]. 

Then we show that the stretched exponential can be considered as a fairly good 
approximation for stress relaxation in a composite with delayed power law and non-zero 
impact response to a stress step, after an initial time interval and during a significant 
part (but not towards the end) of stress relaxation. We derive the same value for the 
parameter of the hyperbolic response and for the parameter of the stretched exponential, 
as well defined function of the basic parameters of the viscoelastic hierarchy. 

To overcome certain limitations of the functional equation approach, we construct a 
fractal transmission line model for the viscoelastic hierarchy of the composite solid. We 
suggest further analysis of the strain and stress dynamics using this tool, originated in 
the electrical engineering field and successfully applied to the analysis of constant phase 
angle impedances of electrode interfaces and low frequency fractional power responses 
in dielectrics, by Le Méhauté, Nigmatullin and Nivanen [13]. 


MATHEMATICAL MODELING 


' A preliminary and partial version of this approach, in the framework of the mathematical 
modelling of mechanical vibrations of structures and machines, under the title “Vibrations of 
viscoelastic, power-law solids”, was presented by one the authors in the First Regional Seminar 
of Latin America and Caribbean on Vibration Analysis and Diagnostics, Bariloche’s Atomic 
Centre, Bariloche, Argentina, September 14 to 18, 1987. 


Let us consider a biphasic composite material. One phase is a linear viscous matrix. The 
other phase is a collection of linear elastic aggregates composed by many elements of 
different sizes and shapes. 


We apply a normal stress o,=o0n the opposite faces of a volume element. Let 
€, = €be the resulting strain in the direction of the applied stress. (Figure 1) 


Figure 1.Element of volume of a biphasic composite with hierarchical structure. 


Working in the framework of the elementary theory of composite materials, we suppose 
that the connection between stress and strain can be obtained in one dimension, for this 


composite, using a sequence of nested viscoelastic boxes (Figure 2). 
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Figure 2. Viscoelastic behaviour of the composite represented by a set of nested viscoelastic boxes. 


The viscoelastic box of order n is composed of N,+lelements in parallel: one linear 
viscous element with damping parameter 7,and N, complex viscoelastic elements. 
Each of these complex elements is the in series combination of a linear elastic element 
with compliance J, and a viscoelastic box of order n+1. 

Then, ifo,, (t)is the stress applied to the viscoelastic box of order n at time?, and if €,, (t) 
is its corresponding strain response, we have: 


,(t)= [ 0,(e-1')-0,(¢)-de (1) 


We assume that the stress and strain are always zero when t<0. The creep kernel of the 
viscoelastic boxes of the level n in the hierarchical viscoelastic structure is Q, (t). 


The viscoelastic box of order 1, which contains all the other boxes nested in its interior, 
is supposed to be in series with a single elastic element of compliance J,. This element 
gives the impact or immediate elastic response of the considered element of volume of 
the composite material. If ¢', (t) is the strain produced in the elastic element of 


compliance J,,, then: é,(t)=e", (t)+é,,,(t) (2) 


But é', (t) = J, ; On (1) » SO that: é,, (1) <i J, ° On (r) zg [ Ora (t ae ') : On (r') : dt' (3) 


But in parallel to the combination of elastic compliance with viscoelastic box of order n 


0 
ao (t) (4) 


If g,,,is the volume fraction of this viscous element, referred to the whole volume of box 


there is a viscous element. The stress acting through it verifies: o,,,(t)=~7, - 


number n, and if @,is the volume fraction of the in series combination of an elastic 


element with box number n+1, from the elementary theory of composite materials we 
have, applying the popular rule of mixtures: 


o,(t)=N, “P, Onalt)+Q,y -,,(t) (5) 
The volume fractions verify, for every n: N,,°Q,+Qy =1 (6) 


The possibility of applying the rule of mixtures, at least as a convenient approximation, 
can be justified as follows. If A, is a cross sectional area corresponding to the boxes of 


level n in the viscoelastic hierarchy and A,,is the area corresponding to the dashpot, 


then: A, -o,,(t)= N° Ana? S, 


n+l n+l 


n 


(+A, -0,.(0) 80 0,()=N, Ae, 4(0)+ ol 


n n 


AsA,=N,-A,,,+A,,, 80 we have ny, Anat 4. So oy 


n+l nv? A, A, 
However, as was discovered by the French geologist Achilles Delesse already in 1842, 
the average fraction of the area occupied by a certain component in a set of parallel 
cross sections of a composite body can be used as an estimation of the volume fraction 


of the component in the body [14]. As consequence of this connection between area 


: ‘ : pea tear le) 
averages and volumes of a given phase in a body of composite material, if —** is 
interpreted as an estimation of the average fraction of parallel cross sectional areas of 
viscoelastic boxes of level n+1 inside a box of level, taking into account Delesse’s rule, 
A 


nv 


A Nectse : 
we have 9, = a An analogous derivation gives @,,, = = 


n n 


Of course both stresses 


and strains must be interpreted also as suitable averages. 


Derivation of a functional equation for the creep kernel Q,(r) 
Taking Laplace’s transforms, and defining f(s)= [- e". f(t)-dt, it follows from (1), 


(3), (4) and (5): &(s)=0,(s)-6,(s) (7) 
G,,(s)= 5-7, -,(s) (8) 


é,(s)=J,,° GaAs )+0,.(s ): G,,,(s) (9) 
G,, s)= No "PD, ° 9, Culs )+9,,° G,,(s ) (10) 
Eliminating g,(s), &,,,(s)andG&,,(s) between (7), (8), (9) and (10) we obtain: 
= 1 
= 11 
Q,(s) > pn Nee. (11) 
Be Os) 


Now, if we assume an infinite hierarchy of nested viscoelastic boxes, then from (11) we 
obtain the following continued fraction: 


— 1 
Q(s)= 7 (12) 
PA 
Py Me S+ 
1 
dit e 
2'Pr 
MS + 
Pr, °My ae 


Now, let us suppose that: VN, =N J,=J 9,=9@ @Q,, =@,, for every n. 
Then (6) reduces to N-9+@,=1 (13) 


If the viscous dissipation of elastic energy is produced mainly in the region of matrix 
surrounding each aggregate-matrix interface, we can expect that as the aggregate 
becomes smaller, the viscosity becomes larger (the other conditions remaining 
unchanged). So, introducing a parameter £ between O and 1, we assume that the 
viscosity coefficients of the different levels of the viscoelastic hierarchy verify: 


n—| 
m (oS =m, (14) 
So, as the volume fraction g is less than 1, 7, grows scaled by = from each level to 
? 


the following level of the hierarchy. Then (12) reduces to: 


Z 1 
O,(s)= - (15) 
2 

Py, Mm S+ 1 
J+ N. 
ples f 
? J+ 
pel eet 22 
v 2 
(0°) Jt... 


Taking (14) into account, and splitting the volume fraction g = vy’ -g"”,, it is possible to 
re-cast the continued fraction (12) as follows: 


= 1 
O,(s)= N-o? (16) 
PM S+ J i 1 
gp” ee N-g'? 
v 1 J 1 
pp 
(0 } PN St... 


As shown in the mathematical Appendix, this expression is convergent because NV -g— <1 


For a given set of parameters of the viscoelastic hierarchy, (16) defines a function of s 
for positive values of the Laplace’s variable. 


If we substitute in (16) s by g’-s we obtain, after some operations, the following 
functional equation: 


0 (a? 
J 
oth 
1+ OKs) 
mi 
The characteristic time T is defined by the equation T=9,-7,-:J (18) 


Equation (17) can be solved numerically. However, here we want to proceed analytically 
as far as possible, so we are going to derive an asymptotic functional equation whose 
analytic solution is known. As also shown in the mathematical Appendix, assuming that 


n-g'? >1, when s\O then Q(s)Too and s-O,(s))0 As consequence, when Ts is 


small enough relative to 1 we will have 2a) small relative to 1 and ats large 
relative to 1, so (16) may be approximated by: 

= —— 

B 
-sJ=—— -OQ|s 19 

Ole’ -s)=— Bs) (19) 

As shown in the mathematical Appendix, the regular solutions of this functional equation 
~ C 
are O(s)= a (20) 
The parameter p is given as the following function of some parameters of the viscoelastic 
5 ) 

hierarchy: ee (21) 


ay 


As will be shown later, p is a parameter comprised always between 0 and 1. 


Inverting the Laplace transform, we obtain the desired delayed strain kernel that gives a 
power law strain response to a step of stress input: 


t a 
a,()= 2,04) 2) (22) 
Here a@=1-—p (also comprised between 0 and 1), and ¢, is an arbitrary but positive 
reference time. If we taket, =T, we can expect that (21) approximates the viscoelastic 


strain kernel fort >7. However, C remains undetermined, as well as Q,(t,) which is 


proportional toC. To overcome this limitation without recurring to (17), let us try a 
different approach to the problem. 


Derivation of a transmission line equation for the viscoelastic hierarchy 
From (5) we have: 2, °C. (t) 7 o,,(t)— N° Pn? Gna (t) (23) 


nv 


From (23), taking into account equations (2) and (4), that is é,(t)=e", (t)+ é,,,(t) 


0 eae ’ 
ando,,, (t)= 71, SF eult), as well as the constitutive relation between stress and strain for 
t 


the in series spring inside boxes of level n, that is’, (t)=J,-o 


n+l 


(t) , we obtain the 
following system of linear differential equations: 


8. (p) = Enalt)—4.(0) 

Prv My = Enlt) N,Q, 
The idea now is to approximate (24) by a partial differential equation of the transmission 
line type used in electric engineering, define suitable boundary conditions and obtain a 
solution that asymptotically gives the power law for strain but without undefined 
multiplicative constants. 
Now, we take the Laplace transform in both members of (24): 
(é,.1 (s)-, (s)) N,-@,: (é, (s)- Ent (s)) (25) 


J n n J 


n-1 n 


After some reordering, (25) may be put as follows, beingT, = @,,,-77, -J 


Kale) teal) (24) 


n-l n 


Pry An -§-&,(s)= 
nds 


1B )=B, (0)-2-5(9)4 (so Me BZ) (6, (9)-F(6) 


n 


If the branching numbers, volume fractions and elastic compliances are constants then: 
T, +8, (s)=,4(s)—2-2,(s)+&,4(s)+(I-N-@)-(2,(8)-&,.(s)) (26) 
Following the approach of Le Méauthé and co-authors in Chapter 7 of their already cited 


book, let us define ¢ = ms and introduce a new variable: x= sae a (7) 
N:-o N-o 


1 n 
If the parameters of the dashpots grow like77, = (5) -77, taking into account that, from 
Pp 


4 nin ae -Inx 
(27), n= + we obtain: (5) =e | aig VA esta (28) 
QP 
-p 
Here, by definition: pe nl -o'*) (29) 


As we shall see in the discussion, v is always positive. 
If we take nas a continuous real variable, then x = €” is a regular real function of7 . 


We define: &, (s)=2(x,s)=&(E",s) Then £.s()-8,(9)=3{ 2.5} ACs) may be 


approximated by ad 2 a(3.s] é,(s)—é,,,(s)=&(x,s)-&(E-x,5) may be 


approximated by—.x-(é—1): < &(x,s) and &, ,(s)—2-é,(s)+é,,,(s) may be 
x 


approximated by x-(é ys +2 E(x, | 
x x 


Taking all this into account, after several operations we obtain the following transmission 
line equation in Laplace’s domain: 


Sas) PPE AU (30) 
Ox 
1 n,-J-N? 
Here, with € = Veo ; ee ete (31) 


Equation (30) allows an approximate description of a finite viscoelastic hierarchy, say of 
N levels. This differential equation is equivalent to a modified Bessel equation and from 
a mathematical point of view is identical with the self-similar transmission line equation 
derived by Le Meauthé and coauthors [13]. 

The solutions of equation (30) are linear combinations of modified Bessel functions, 


ay 
both of the first kind x -J {2-pas and of the second kind 


1 1 In 
2 Dp et 1 [to 
x Kear i: The parameter p verifies: p= = 


ao (32) 
+y (2 


gp’ 


So, it takes the same value given by equation (21) and obtained from the solution of the 
approximate functional equation (19). 

When 7,-s is positive and much less than 1, the dominant term in the formula of the 
Laplace’s transform O,(s) for the viscoelastic creep kernel Q, G ) can be estimated by a 
suitable modification of the corresponding formula for the input impedance of an infinite 
self-similar transmission line: 


O(s)=J+p'??- ey & y (33) 


Here, for x > 0, T(x) = | e“-u*'-du is the Gamma Function. 
0 
As consequence, the constant C in formula (20) is now determined and can be written 


= pn 2P., T(p) P 1 : 
eas (a =] ae 


€ 


Derivation of the stretched exponential stress relaxation 
If a linear viscoelastic material has both instantaneous (elastic) ¢,(t) and delayed 


a (t) strain response to an applied stress o(t)that is zero for negative times, the total 


strain &,(r)+é,(t) is given by: é(t) = J, -o(t)+ [ K(t-1')-o(t')-dt' (35) 


If the case of the composite material studied in this work, the creep kernel K (t) may be 
identified with Q, (t), the kernel of the box in level one in the viscoelastic hierarchy. 
Inverting (35) and introducing an elastic impact modulus VM, such thatJ,-M,=1 anda 
stress relaxation kernel L(t): o(t)=M, -e(t)- i L(t -1')-e(t')-dt' (36) 
Taking Laplace’s transforms in (35) and (36), it is possible to relate the transform L (s) 


of the stress relaxation kernel with the transform K(s) of the creep kernel: 


(37) 


C 


Now we substitute the Laplace’s transform O( a a of the power law kernel (see 


Ss 
: ; > Cc 
equation (20)) in place of K (s ) in (37), rs [z I and we obtain: 
T., 


ny, (a) ae 


1 1 

Dp ae a) = 

From x = 2) and from (34) it follows that V= cal -p? { r(p) i (39) 
0 T J rl - p) 


Inverting the transform (37): Lt) _ x { fe sy) (1) (2 | (40) 
T 


Mo. fe rtqXt “41 ((n+1)- p) 


As usual, T(z) represent the Gamma function. An equivalent expression was called 


7 € 


fractional exponential function by Rabotnovy, already in 1948, since for p=1 reduces to 
t 


£}.6(%) 


the decaying ordinary exponential. In our case, when p =1 we obtain: 


If, as mentioned in [15], the function ered is approximated by 
mo 1 (\n 


AF 
From (36) the relaxation of stress at a constant strain €, for a viscoelastic body is given 


€ 


es —p(2-p)’- at 
p-(2-py’ eR pyu? ‘hen: L(t) M, 2 {4e) p22 vy is p(2-p (24) (a1) 


by o(t)=M(t)-€ with, by definition: M(t)= (™, - [ L{t’'): ar) (42) 


Taking into account (41), from (40) it follows that the stress relaxes according to a 


LY 
stretched exponential: o(t) =M (t): Ey =O" te (43) 
In this equationo, = M,-& 

From equations (41) and (39) we obtain the characteristic stress relaxation time T, 
given as function of the characteristic time 7, of the hyperbolic strain response: 


t, = & (44) 


ema er( tO)! (oF 


(I— p J 


The asymptotic percussion stress relaxation mode in the viscoelastic hierarchy 
The stretched exponential is exact forp=1. Considered as an approximation it 


progressively weakens when p decreases towards zero. However it is a fairly good 


approximation for p > Vy, when Ve is greater than one but not too much. When 


vA — +00 it is possible to derive the following asymptotic expansion” for the stress 


relaxation in response to a step in strain [14]: o(t) Fo yi ‘Ngee T(np)-sin(z-n- p) 


iS x . np 
Tt, 


So, for ue big enough, o(t) x Fo. T(p)-sin(z-p) (45) 
: ely) 
wh, 


DISCUSION 


The parameter p of the power law strain response 


and es are 
? 

greater than 1. But from both, experimental results and thermodynamic requirements 

related with elastic energy dissipation, we know that p must be less than 1. This is 


From (32) we conclude that p must be positive, because both € = 


equivalent to the restrictionv>1. Taking into account the definition of v given in 
equation (29), as well as the inequality N-@ <1, we see thatv >1 is verified if and only 


if N-g'’ >1. The physical origin of this last inequality and the physical origin of the 


scale law for the dashpot viscous moduli 77, (9° yn =7,, can be found examining the 
geometric scaling laws of the viscoelastic hierarchy. 

Let us consider an element of volume V, in the composite material corresponding to a 
combination in series of an elastic element and a whole viscoelastic box of level n. 


* It diverges for every time instant. 


Let us introduce a characteristic length/, corresponding to volume V, and a 


n 


characteristic area S, ~/? of the interfaces between the volume element and its viscous 


| 2 
matrix environment. We assume that/,,, = oy ‘1, so S44, & gy -S,, and V, 


n+l P VY, 
Furthermore, we suppose that the force F taken by the viscous matrix in the volume 
OE, 


Ot 


element is given by F, =w-S, - Then, the force per unit volume o,,, is given by: 


We suppose now that the viscous stress ois related with the density of volumetric 


= 0€ 
viscous forces by a constant factor: o,,=y-:o0,, So thato,, =7,,- a” with 
t 


1 n-l 
at S 
n=[0" | ‘7, and mee al waar 


1 


i 7 a te : n—| 
As consequence, we obtain the scale law of viscous dissipation 77, (9? y = 7, with 


p => According to this model, the parameter of the power law and of the stretched 


1 
or ) 
exponential verifies p =————_<.._ Now v is positive and p is less than one if 


of ° 


ee = N 0” = (w pg? ) is greater than |: this happens if the total area related 


with viscous dissipation increases from one level of the hierarchy to the next one. 


The relaxation times in a self-similar viscoelastic hierarchy 

The self-similar nature in physical space of the viscoelastic hierarchy proposed here 
gives rise to power law responses, and power law responses point to scale invariance in 
time. Scale invariance, at its turn, implies a manifold of time scales in the delayed strain 


leu (t) 


= J(1)-t? -o,, such that the local ti I ia (al 
responses &,(t)=J(1)-t? «oy, such that the local time scale lde,,()/ae 


let) t a Mente 
——~~___ = —. The characteristic times 7, 
de, (t)/ar| P 


of the response 


at the instanttis always proportional tor: 


do not appear in the local time scales. 


In the case of stretched exponential stress relaxation ¢ (r)~ o(0)-exp |- [ C/A 
r T, 
o,(t) 


ldo,(t)/dt] 
local time scale of the delayed strain response. But now the characteristic time 


y] the local 


Pp 
time scale is: "6 -t It grows liket'”, slower than the growth of the 
t g g 


T, appears in the formula for the local time scale: the stretched exponential does not 
have scale invariance in time. 


CONCLUSIONS 


(a) A unified phenomenological approach to power law strain and _ stretched 
exponential stress relaxation was developed in this paper, for a simplified 
example of viscoelastic hierarchy. 

(b) An analytical formula for the power law parameter as function of the branching 
number and volume fractions of the hierarchical composite material was derived 
for the simplified example. In this case the parameter of the stretched 
exponential turned to be the same as the power law parameter. 

(c) Formulae were obtained also that allows us to calculate the characteristic times 
t,and t, as well defined functions of the mechanical parameters of the 


composite material. 
(d) Due to its scale invariance, the hyperbolic strain response has an infinite and 


: t tek, ate 
unbounded set of local time scales—. The characteristic time t, must not be 
p 


considered as a time scale of the hyperbolic response. It gives a measure of the 
numerical order of the interval of time (since the application of a step in stress) 
after which the power law can be applied. 

(e) The characteristic time 7, gives a measure of the numerical order of the interval 


of time (since the application of a step in strain) after which the stretched 
exponential law for stress relaxation can be applied. The stretched exponential 


P 
also has an infinite and unbounded set of local time scales ("% ) -t Contrary 


to the hyperbolic response case, now 7, appears in formula for the local time 


scale of stress relaxation. 

(f) If properly validated in practice, the kind of approach developed here could be 
used as a guide to the design of power law viscoelastic materials with desired 
values of the power law parameter and characteristic timest, andT,. 


(g) In principle the equation (12) for the Laplace’s transform of the creep kernel 
allows to take into account random branching numbers as well as random elastic 
and viscous moduli. The resulting equation could be approximately solved 
applying a method developed by Kaplan and Gray to study the effect of disorder 
on a fractal model for the alternating current response of rough interfaces [16]. 

(h) The transmission line approach enables us to model different in series and in 
parallel viscoelastic hierarchies, considered as generalized quadrupoles. In this 
case strain corresponds to voltage and mechanical force (stress multiplied by 
area) corresponds to current. Dashpots correspond to capacitors and elastic 
elements correspond to resistances. With these correspondences, the well- 
developed mathematical tools of distributed parameter electric networks can be 
applied to study viscoelastic hierarchies. 

(1) Finite viscoelastic hierarchies show, besides the high frequency cut off that 
corresponds to the beginning of the deformation process, a low frequency cut off 


for long times of deformation. In the model based in a set of nested boxes, the 
low frequency cut off behaviour depends of the termination (elastic elements or 
dashpots) of the last box in the hierarchy. 


MATHEMATICAL APPENDIX 


1. A continuous fraction like Q,(s) originates a sequence of the so called reduced 


fractions: 
1 
= 1 1 
R)=——— R,)= az Rls)- — 
?,° m° oS cts Pe 2, 7, ar oo 1 
g” eg” em 
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1 
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When s is positive: 0 < R,(s)< R,(s)< R,(s)<...< O,(s)<...< R;(s)< R,(s)< R,(s) (A2) 
If R,,,(s )- R (s) — 0 whenn — +oofor every positive s , the continuous fraction defines 


a function of s that is a solution of the functional equation (17). See [17], Chapter 2. 
1 1 


(T : s)* : (v | ot 


is asymptotically dominant, Ry als )- R,,(s) tends to 


In our case, withT =@,-77,-J wehave: 0< Ry,.,(8)— Roy (s) < 


1 
Because N -g < 1 and ————_, i 
(v-o)""” 
zero whenk —-+00. As consequence the functional equation for the Laplace’s 
transform of the creep kernel has at least one positive solution. 

2. Let us consider this positive solution Q,(s)of the above mentioned functional 
1 N-o'? 


equation, rewritten: = =T-s+ (A3) 
gp Oo? -s)/J OG ys 


From (A2): R,(s)<@,(s) From (A1): R,(s)> >0 whens +0. So, the lower 
oe, 


limit Q, of Q,(s) when s\ 0 must be positive or +o. 
Furthermore, liminf O, (9° . s) =Q, 
From (A3) and taking into account the properties of upper and lower limits (see [18]): 


1 _ N-@ (A4) 
Q,/J  14+Q,/J 
Because N -g <1 the only admissible solution of (A4) is: QO, = +00 (A5) 
Butif liminf,, Q,(s) =+00 then: lim OQ, (s) = +00 (A6) 


3. Now we consider the functional equation rewritten again thus: 


s-g° O(o*-s)_ 1 (A7) 


s-Q(s) 


S+ 


J 
From (A7) it follows that either @, =limsup,,, s-O,(s)=0 or 4, is finite and verifies: 


AW spr? 
£.= oo Seow Then: £,= ees But if N-g'’ >1 this last equation 
N-g'? T 
T+ 
ee 
implies that ? must be negative, which is not possible. So, the only possible solution is: 
¢, =limsup |, s- O(s)=0 As consequence: lim 155° O(s)=0 (A8) 
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